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ABSTRACT 

There are only two methods for estimating the mass distribution in the outer regions of galaxy clus- 
ters, where virial equilibrium does not hold: weak gravitational lensing and identification of caustics 
in redshift space. For the first time, we apply both methods to three clusters: A2390, MS1358 and 
CI 0024. The two measures are in remarkably good agreement out to ~ 2/i _1 Mpc from the cluster 
centers. This result demonstrates that the caustic technique is a valuable complement to weak lensing. 
With a few tens of redshifts per (/i _1 Mpc) 2 within the cluster, the caustic method is applicable for 
any z < 0.5. 

Subject headings: cosmology: miscellaneous - cosmology: observations - galaxies: clusters: individual 
(Abell 2390, CI 0024+1654, EMSS 1358+6254) - gravitational lensing 



1. INTRODUCTION 

The relative distributions of mass and light in the uni- 
verse have remained a profound and central mystery in 
cosmology for more than seventy years. Since Zwicky's 
pioneering use of the virial theorem to di scover dark mat- 
ter in the Coma cluster l)Zwickvlll933t) . the range and 
sophistication of methods for estimating cluster masses 
and mass profiles have increased to include a host of dy- 
namical measures, X-ray estimates and strong and weak 
gravitational lensing determinations. 

Different mass estimators applied to rich clusters of 
galaxies constrain the mass distribution on different 
scales. Strong lensing generally provides constraints on 
very small scales (< O.l/i -1 Mpc). Virial mass esti- 
mates, including Jeans' analysis, assume dynamical equi- 
librium and apply only within the virial radius. Mass 
estimates based on X-ray observations assume hydro- 
static equilibrium a nd rarely extend beyond one-half 
of the virial radius j^^ye rowicz. Neumann, fc Reinricbl 
12002 IPratt fc Arnaudll2002f) . 

At larger clustrocentric radii where equilib- 
rium assumptions break down, there exist only 
two techniques for mass estimatio n: weak lensing 
ijKaiser. Squires, fe Broad hurstj|1 995[) and the redshift- 
space caustic technique ( Diaferio & Gclle rll9'9'^: iDiaferiol 
119991 D99 hereafter). Both techniques enable determi- 
nation of the mass distribution from the cluster center 
to distances larger than the virial radius. 

The ca ustic technique ha s been applied to many local 
clusters ijRines et alJ 120031 and references therein). At 
small clustrocentric radii, caustic estimates agree well 
with the traditional virial analyses in the optical and X- 
ray bands. At larger radii, the caustic technique is still 
valid, but its mass estimates were tested against TV-body 
simulations only (D99). 

Here we discuss the first comparison of mass estimates 
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from the caustic technique and weak lensing. Only re- 
cently have sufficient lensing and spectroscopic data be- 
come available to make this comparison. Both techniques 
have known systematic uncertainties: these comparisons 
test the importance of these systematics. 

In this Letter, we examine mass profile measure- 
ments for three intermediate redshift clusters: A2390, 
MS1358+6254, and CI 0024+1654. 

2. THE CAUSTIC TECHNIQUE 

Cluster galaxies plotted in a redshift space diagram 
(line-of-sight velocity v vs. projected distances R from 
the cluster center) distribute in a characteristic trumpet 
shape. The boundaries of this trum pet are called caustics 
l)Kaiserll987tlRegos fc Gellerlll989ft . By assuming spher- 
ical symmetry and hierarchical clustering for the forma- 
tion of the large-scale structure, the caustic mass esti- 
mator relates the caustic amplitude, the trumpet width 
in v at each radius R, A{R), to the escape velocity from 
the gravitational potential well generated by the cluster. 

The procedure developed by D99 provides an auto- 
matic method for locating the caustics and determining 
their amplitude. First, the procedure arranges all the 
galaxies in the field in a binary tree and finds the clus- 
ter members. The cluster members determine the cen- 
ter of the cluster, its one-dimensional velocity dispersion 
(v 2 ) 1 / 2 , and its radius {R), the mean projected distance 
of the members from the cluster center. Table H lists 
these quantities for the three clusters. 

The procedure next determines the threshold k which 
enters the caustic equation f q (R,v) = k. Here, f q (R,v) 
is the galaxy density distribution in the redshift dia- 
gram, smoothed with an adaptive kernel. The param- 
eter q sets the scaling between the quantities R and v. 
We choose the parameter k by minimizing the quantity 
S(k,(R)) = \(v* sc ) Kj{R) -A{v 2 )\ 2 , where {v 2 csc ) K , {R) = 

f( R) A 2 (R)ip(R)dR/ f( R) cp(R)dR is the mean caustic 
amplitude within (R) and f(R) — J f q {R,v)dv. 

D99 shows that the three-dimensional cumulative mass 
profile can now be estimated as 

1 r 

GM(< r) = - A 2 {R)dR . (1) 

2 Jo 

The error bars on individual data points are propor- 
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TABLE 1 
Cluster parameters 



cluster 


FOV (a x S) 


Afield 6 [zi, Z 2 ] 


N 


Q 


5 


z 




R 


r 3 


c 


^200 


A2390 
MS1358 
CI 0024 


43f8 x 7.'4 
20.'9 x 21.'4 
20.'0 x 24f3 


351 G [0.1,0.4] 
360 6 [0.1,0.5] 
399 e [0.3,0.5] 


210 
282 
251 


21 h 53 m 35!53 
13 h 59 m 49!65 
h 26 m 35?90 


17° 42' 3'.' 16 
62° 30' 55'/ 87 
17° 9'41?10 


0.2284 
0.3289 
0.3941 


1154 
996 
937 


0.85 
0.76 
0.74 


0.14 ±0.17 
0.14 ±0.09 
0.12 ±0.11 


11 ± 12 
7.7 ±4.3 
8.6 ±7.7 


1.5 ±2.4 
1.1 ±0.9 
1.0 ± 1.3 



Note. — Col. 2: Field of view (FOV); Col. 3: no. of galaxies in the FOV within the redshift range [21,22]; Col. 4: no. of members; Cols. 
5-7: cluster center coordinates; Col. 8: cluster velocity dispersion; Col. 9: cluster size; Cols. 10-12: NFW fit parameters. Celestial coordinates are 
J2000. Velocities are in km s~ , lengths in h^ 1 Mpc. 



tional to the inverse of the galaxy number density within 
the caustics (D99). This recipe quantifies the uncer- 
tainty in the mass estimate which mostly results from 
deviations from spherical symm etry. The recipe was cal- 
ibrated on TV-body simulations l)Kauffmann et af]ll999l) 
that generally showed less cleanly defined caustics than 
in the real Universe. Therefore, we suspect that these un- 
certainties are smaller for real clusters than in the simula- 
tions. The small scatter (< 30%) around the equivalence 
relati on between X-ray and caustic masses ijRines et al.l 
2003) suggests that the simulations indeed overestimate 
the errors in the caustic technique at small radii. If 30% 
represents a rough estimate of the correct caustic mass 
uncertainty at all radii, the D99 recipe typically overesti- 
mates this uncertainty by a factor of two. Nevertheless, 
because it is the only available prescription for evaluat- 
ing the error, we use the conservative D99 prescription. 
Comparison of gravitational lensing and caustic measure- 
ments for large samples of clusters in the redhsift range 
0.2-0.8 will test the accuracy of this recipe. 

3. MASS COMPARISON 

Mass profile estimates of high-redshift clusters de- 
pend on the assumed cosmological parameters: phys- 
ical distances, X-ray and weak lensing cumulative 
mass profiles scale as the angular diameter dis- 
tance Da- Moreover, i f one derives a best-fitting 
iNavarro. Frenk. fe White! l)1997|) (NFW) density pro- 
file p(r,z) = 5 c p C at{z){r/r s )- l (l + r/r s )~ 2 , with 
pait(z) — 3H 2 (z)/8itG the critical density of the Uni- 
verse, h 2 (z) = H§[n {i+z) 3 +(i-n -n A ){i+z) 2 +n A }, 

S c = c 3 (200/3)[ln(l + c)-c/(l + c)]- 1 and c = r 200 /r s the 
concentration parameter, c also depends (non- linearly) 
on Da, because 5 c p CT it(z) scales as D^ 2 . Below, all 
quantities assume fio = 0.3, Q\ = 0.7 and Hq = 
100/1 km s" 1 Mpc" 1 . 

Figure shows the redshift diagrams of the three clus- 
ters with the caustic location (upper panels) and the 
mass profiles estimated with the caustic technique, grav- 
itational lensing and X-ray data (middle and lower pan- 
els). Gravitational lensing measures all the mass pro- 
jected onto the sky along the line of sight. Therefore, 
we distinguish between three-dimensional (middle pan- 
els) and projected (lower panels) cumulative mass pro- 
files. Radial distances are three-dimensional (r) or pro- 
jected onto the sky (R). 

The solid lines in Figure Q show the best-fitting NFW 
profile with parameters listed in Table ^ To com- 
pute these fits, we only used the data points within 
Him = l/i -1 Mpc, a conservative radius beyond which 
the NFW mass profile might not be a good description 



of the actual profile. For all clusters, the data points be- 
yond 1/t, -1 Mpc do indeed agree with the NFW model, 
indicating that the correct choice of ru m is irrelevant. In 
any case, the fit parameters and their errors are only 
indicative, because the individual data points are corre- 
lated. Moreover, the NFW fit parameters are correlated 
even with independent data points. Keeping one of the 
two parameters, c or r s , fixed in our fits reduces their 
relative errors to ~ 10%. 

For each cluster, we also show the best fits determined 
from the weak lensing (dashed lines) and X-ray (dotted 
lines) measurements. We now comment on each cluster 
separately. 

A23 90 is a rich cluster at z = 0.228 w ith op- 
tical llLe Borgne et all I199H lYee et alJ 119961) X-ray 
l|Bohringer et alJ Il998t lAllen. Ettpri. fc Fabianl 120011) 
and both weak JSauires et alJ 1996|) and strong 
jPello et allll991i iPierre et alJl!996j) gravitational lens- 
i ng observation s . 

Squi res et al.l l|1996l) compare the weak lensing data 
within ~ 260" with a singular isothermal model with 
velocity dispersion a — 1093 km s _1 taken from 
iCarlberg et alJ l)1996j) . The isothermal model underpre- 
dicts the amount of mass actually measured in the range 
0.46 — 0.67/i _1 Mpc (left-bottom panel in Figure^); how- 
ever, this model is in good agreem ent with th e best- 
fitting NFW mass profile derived bv lAllen etaD (2001) 
from Chandra observations. They find r s = OAA^J^h' 1 
Mpc, c = 3.6±?;g and r 2Q0 = 1.6±?;?fa- 1 Mpc. 

By using the galaxy redshift survey bv lYee et alJlTl 996) 
and by assuming dynamical equilibrium. iCarlberget alJ 
(1996) estimate a mass M(< 3.3/i _:L Mpc) = (2.7±0.4) x 
10 15 h- 1 M e . The caustic mass (1.4 ± 1.2) x lO 15 ^" 1 !^ 
and the mass 1.8 x 10 15 /i _1 Mq extrapolated from the 
weak lensing isothermal model are 48% and 33% smaller 
than this virial mass, but within its 3-cr uncertainty. 

At s maller radii. A23 90 sports spectacular arcs and 
arclets llPello et alJll991lh some of which have measured 
redshi ftsllBezecourt fc SoucaiHll997t iFrve fc Broadhurstl 
Il998t iPello et alJ 119991) . IPierre et alJ l|1996ft use the 
brightest strongly lensed arc and its surrounding shear 
to derive the projected total enclosed mass M(< 
97/i~ 1 kpc) = (8.0 ± 1.0) xlO^/i^Mo 4 (solid diamond 
in Figure Q), in agreement with the mass (1.2 ± 0.7) x 
1O 14 /i _1 M implied by the projection of the NFW fit 

4 In this Letter, we rescale each strong lensing mass found in the 
literature by the effective lensing distance D[Di s /D s appropriate 
to a universe with f2rj = 0.3 and Q\ = 0.7; D;, D a and D; s are the 
angular distances to the cluster, to the source of the lensed image 
and between the cluster and the source, respectively. 
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Fig. 1. — The left, middle and right columns are for A2390, MS1358 and CI 0024, respectively. Top panels: Redshift diagrams with the 
galaxies (dots) and caustic locations (solid lines). Line-of-sight velocities v are in the cluster rest-frame. Middle panels: Three-dimensional 
cumulative mass profiles. The solid squares show the caustic mass estimates; the solid lines are the best-fitting NFW profile s to the data 
points within lh —1 Mpc; the dotted lines are the best-fitting NFW profiles to the X-ray me asures (from left to right: IAller^^alJ|200 
Arabadiis ct al. 2002; IQta et aH2004l): the da shed lines are the best-fitting isothermal (A2390, Squires et al. 1996; MS1358, Hoekstra ct al 
11998) or NFW models fCl 0024. IKneib et alJI2003n to the gravitational lensing measures. The left and right vertical dotte d lines show the 
radius of the X-ray and gravitational lensing fields of view, respectively. The two filled circles show the virial estimates bv lCarlberg et all 
(1996) of A2390 and MS1358. Bottom p anels: Projected cum ulative mass profiles; lines are as in the middle panels. The open diamonds 
show the weak lensing measures: A2390, ISauires et alJ 119961): MS1358, lower limit bvlHoek stra et alJ 119981) to the mass p rofile. Filled 
diamonds show the strong lensing measures: A23 90, IPierre et alJ 11999): MS1358: l Allenl 119981) from the measurement by IFranx et all 
11997); CI 0024: upper symbol. ITvson et alJ 119981) . lower symbol. iBro^dhurst et al.l 12000?) . Error bars in all panels are 1-cr; error bars on 
points where they seem to be missing are smaller than the symbol size. 



to the caustic mass; the strong lensing mass also agrees 
with the mass 8.5 x 1O 13 /i _1 M implied by the weak 
lensing isothermal model, and is just above the 68% con- 
fi dence bound deriv ed with the X-ray analysis (Figure 8 
o flAllen et aLpOOl 

IPierre et al.l l|1996[) derive the strong lensing mass by 
assumi ng that the arc is a single lensed galaxy at z = 
0.913. iFrve fe Broadhurstj l)1998j) later showed that the 
fainter part of this arc actually is a second lensed galaxy 
at z = 1.033. The redshifts of the arcs and arclets, which 
are available now but not at the time of Pierre et al.'s 
analysis, urges a reformulation of the lensing model of 
the core of A2390. However, we expect that a newly de- 
rived mass will no t substantially differ from the mass of 
IPierre et alJ l|1996[K because the mass estimated with the 
simplest lensing models, which provide the most inaccu- 
rate measures, probably are within 30% of th e true value 
ijKochaneck. Schneider fc Wambsgans s 2003). 

MS135 8+6254 is a very rich cluster first discov- 
ered by IZwickv fc HaHogl lfl96l . We collect 381 
re dshifts in the cluster region from the surveys 
of iFabricant. McClintock. fc Bautzl ljl991|) . iFisher et all 
H1998I) and lYee et alJ H19M) ~ 



Hoe kstra et all l)1998|) used HST observations to con- 
struct a weak lensing map of the cluster extending to 
a radius of - 220" = 0.73/1" 1 Mpc. They only de- 
rive a lower limit to the mass profile and find a best- 
fitting singular isothermal model with a = 780 ± 50 
km s -1 (dashed lines in Fi gure H}. More recently, 
lArabadiis. Bautz. fc Garmird l|2002|) analyze a Chandra 
observation of the cluster. They approximate the mass 
profile within ~ 2' = OAh^ 1 Mpc with an NFW pro- 
file, with r s = 88+^Zi" 1 kpc, c = 9.3±'^ and r 20 o = 
Mtnith' 1 Mp c. 

ICaTlbcrg ct al. (1996) assume virial equilibrium to es- 
timate M(< 2.5/i" 1 Mpc) = (1.5 ± 0.2) x I0 l5 h,- 1 M Q 
from their galaxy redshift survey. This mass is more 
than 3-cr above the weak lensing isothermal extrapola- 
tion 7.0 x 10 14 ft- _1 M Q which agrees with the caustic esti- 
mate (6.5 ± 2.8) x 10 u h- 1 M Q . The extrapolation of the 
X-ray fit yields 3.0 x 10 14 /i _1 M Q , a factor of two smaller 
than the caustic mass and a factor of five below the virial 
mass. Probably, the assumption of virial equilibrium at 
this large distance is unrealistic and the extrapolation 
of the X-ray profile, limited to radii < 0.4/i -1 Mpc, is 
unreliable. 
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In the very central reg ion, lAilenlll 1998) uses the strong 
lensing observations by iFranx et alJ (flQQ^h to derive a 
projected mass M(< 69ft _1 kpc) = 4.4x 10 13 ft _1 Af with 
a 20% uncertainty. The projected NFW profile derived 
from the caustics yields a perfectly consistent mass (4.2 ± 
1.3) x lO^ft" 1 M(^. The projected profile s derived by 
iHoekstra et~aTI l|1998j) and lArabadiis et all l)2002fl imply 
the somewhat lower masses 3.0 x 10 13 /i _1 Mq and 2.9 x 
10 13 h^ 1 M Q , respectively. 

The X-ray and weak lensing mass models agree within 
~ 0.8ft -1 Mpc, but underes timate the strong lensing 
mass derived bv lAllenl l)1998[) . The fact that the weak 
lensing mass provides only a lower limit to the mass pro- 
file and the caustic mass is in excellent agreement with 
the strong lensing measurement suggests that the caustic 
mass provides the correct mass profile of MS1358 out to 
- 2ft -1 Mpc. 

Significant tension exists be tween 

lensing iBonnet. Mellier. fe Fort] 119941: 

iTvson. Kochanski. fc deH'Antoniol 11998]) and X-ray 
ijSoucail et alJ 120001: lOta et alJ I2004D mass estimates 
of the CI 0024 cluster. iKneib et alJ $2003(l combine 
their weak lensing measurements from wide field 
imaging with the str ong lensing measurement by 
IBroadhurst et all (|2000|) to derive the best-fitting NFW 
profile with r s = 54 ± 2ft" 1 kpc, c = 18. Ttl'l and 
^200 = 1.0lio"||ft. — 1 Mpc. According to their Figure 12, 
the uncertainty in their mass estimate is always < 10%. 
Our Figure ^ also sh ows the NFW profile which fits 
recent Chandra data l)Ota et al]l2004D . These authors 
derive the NFW profile from a j3- model fit. According 
to its parameters, we find r s = 0.56 ± 0.02ft _1 Mpc, 
c = 1.8 ± 0.3 and r 200 = 1.02 ± 0.18ft -1 Mpc. Our 
caustic estimate lies between the lensing and the X-ray 
fits at r < 0.2/i" 1 Mpc, but it is in excellent agreement 
with the lensing estimate outside ~ 0.5ft -1 Mpc. 

In the cluster central region, there are two strong 
lensing measurements which yield comparable masses. 
However, the very small errors claimed make them in- 
consistent with each othe r: M(< 0.114ft~ 1 Mp c) = 
(1.30 ±0.04) x lO 14 ft- 1 M lIBroadhurst et aljEfflO) . and 
M(< O.ligft^Mp c) = (1.563 ± 0.002) x lQ u h- l M Q 
llTvson et al.l ll998). We scaled the mass reported by 
lTTson"et^uTl^)98fl by assuming z = 1.6 75 for the arc. 
as measured bv IBroadhurst et al J 120001). B v construc- 
tion, the NFW profile of IKneib et alJ l|2003) agrees with 
the former (it yields M(< 0.114ft _1 Mpc) = 1.13 x 
lO 14 ft _1 M0) and therefore disagrees with the latter (it 
yields M(< 0.119ft _1 Mpc) = 1.17 x 10 u ft _1 M Q ). The 
caustic profile gives smaller, but consistent, masses in 
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both cases: (7.9 ± 3.8) x 10 13 /i _1 M© and (8.5 ± 4.0) x 
10 13 ft -1 M©, respectively. The NFW fit to the X-ray data 
yields even smaller masses: 3.8 x IQ^h^ M^ and 4.2 x 
10 13 ft ~ 1 M f n with a ~ 30% typical error. iCzoske et a!] 
(2002) suggest that the peculiar redshift distribution of 
the galaxies within ~ 3.5ft _1 Mpc of the cluster center 
can be explained by a high-speed collision along the line 
of sight between CI 0024 and a less massive cluster. This 
model implies that the X-ray mass estimate based on 
dynamical equilibrium is unreliable. Because the caustic 
and lensing mass estimators are both independent of the 
dynamical state of the cluster, it is reasonable that they 
agree with each other but disagree with the X-ray mass. 
4. CONCLUSION 

For the first time, we compare the only two cluster 
mass estimators that do not rely on the dynamical equi- 
librium of the system: weak gravitational lensing and 
caustics in redshift space. We estimate the caustic mass 
of A2390, MS1358 and CI 0024 within - 2ft -1 Mpc of 
the cluster center. The caustic mass profiles are in very 
good agreement with the lensing profiles. We confirm 
that the discrepancy between lensing and X-ray mass in 
CI 0024 is probably a consequence of the unrelaxed state 
of the cluster which invalidates the X-ray analysis. 

Weak lensing requires accurate photometric wide-field 
surveys in excellent seeing; moreover, the cluster sample 
is somewhat limited to clusters at distances where the 
lensing signal is sufficiently strong. Weak lensing mea- 
sures all of the mass projected along the line of sight, 
resulting in a min imum 20% uncertainty in the cluster 
mass estimates (Ide Putter fc White! 12003) . The caus- 
tic technique, which requires dense wide-field redshift 
surveys, provides a complementary measurement of the 
three-dimensional mass profile of individual clusters at 
moderate redshift; it also yields robust mass profiles for 
clusters in the local universe. 

Future comparison of these techniques for large sam- 
ples of clusters, covering a range of redshifts, will con- 
strain systematic uncertainties in the methods and may 
provide insight into the change in the relative amounts of 
mass in the infall regions and cluster cores as a function 
of lookback time. 
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